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SUMMARY  •  < 

Data  for  Fortran  IV  programs  has  to  be  in  fixed  format.  Some  users  find 
this  constraint  irksome.  Presented'here  a  data  input  scheme  for  Fortran  IV 
which  makes  the  minimum  of  demands  on  the  data  structure.  If  the  data  is  well- 
formed  and  unambiguous  it  will  be  read.  Ill-formed  and  ambiguous  data  will  also 
be  read  and  given  a  reasonable  interpretation,  with  the  location  of  the  suspect 
data  and  the  value  assumed  being  output  as  a  (suppressible)  error  message.  The 
scheme  also  allows  input  variables  to  retain  their  previous  values,  identical 
consecutive  data  to  be  input  in  a  simplified  form  and  alphanumeric  comments  to 
appear  amongst  the  data.  The  whole  data  input  scheme  is  written  in  standard 
Fortran  IV  so  that  the  advantages  of  machine  transferability  of  the  program  are 
retained . 
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]  INTRODUCTION 

Programs  written  in  standard  Fortran  IV  are  limited  to  data  input  with 
fixed  format.  Often  it  is  more  convenient  to  input  data  separated  by  commas  or 
spaces,  for  example,  and  in  recognition  of  this  various  manufacturers  have 
introduced  their  own  'free-format'  data  input  routines.  Unfortunately  these 
are  not  standard  and  programs  and  data  prepared  for  input  in  free-format  on  one 
manufacturer's  machine  cannot,  in  general,  be  transferred  to  those  of  another 
manufacturer.  In  this  Memorandum  free-format  data  input  routines,  written  in 
standard  Fortran  IV  are  presented. 

The  scheme  enables  data  to  be  read  into  various  arrangements  of  variables, 
with  each  routine  starting  to  read  on  a  new  line  of  data  and  continuing  until 
either  a  specified  number  of  numbers  has  been  read  or  an  end  of  read  symbol 
(/  or  $)  is  encountered.  In  specifying  the  routines,  the  aim  has  been  to  give 
the  correct  interpretation  of  all  well-formed  numbers  and,  when  ambiguous  data 
is  encountered,  to  give  a  reasonable  interpretation  and  to  issue  a  warning. 

A  failure  in  the  input  routines  can  only  occur  from  a  system  failure,  for 
example  if  a  number  is  too  large  for  the  computer  to  handle  or  an  attempt  is 
made  to  read  beyond  the  end  of  a  file. 

A  number  of  additional  features  have  been  introduced  for  convenience  in 
using  the  routines.  Alphanumeric  comments  can  be  inserted  into  the  data  by 
enclosing  the  characters  between  inverted  commas  (either  single  or  double) . 

An  input  variable  can  retain  its  previous  values  (that  is,  a  'null'  datum  is 
read)  by  the  use  of  two  consecutive  commas,  optionally  enclosing  blanks.  If 
consecutive  values  being  read  are  the  same,  then  the  data  input  can  be  simplified 
using  the  convention  i *  V  to  denote  i  consecutive  occurrences  of  value  V  . 
The  character  /  ends  the  read  call  forthwith:  any  further  characters  beyond  / 
on  the  same  line  are  ignored. 

Whilst  one  of  the  aims  in  the  specification  of  the  routines  has  been  to 
maintain  compatibility  with  existing  formatted  data  for  Fortran  programs, 
this  has  not  always  been  possible.  For  example,  in  order  to  read  the  two 
numbers  149  and  736  in  Fortran  213  format,  they  would  appear  in  the  data  as 
149736.  The  input  scheme  described  here  would  interpret  this  as  the  single 
number  149736.  For  all  cases  where  formatted  data  is  separated  by  a  non-digit 
character,  the  routines  will  read  in  the  numbers  as  intended.  In  addition,  the 
routines  will  interpret  correctly  free-format  data  prepared  in  accordance  with 
the  proposed  Fortran  77  standard. 
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Details  of  the  data  input  scheme  are  given  in  section  2,  while  the 
programs  themselves  are  described  in  section  3  and  listed  in  Appendix  B.  A 
simplified  users'  guide  to  the  scheme  appears  as  Appendix  C. 

Since  this  data  input  scheme  was  developed,  the  authors'  attention  has 
been  drawn  to  a  similar,  but  less  comprehensive  scheme,  devised  by  D.  Lovell 
which  has  been  in  use  as  part  of  the  data  handling  software  of  the  low-speed 
wind-tunnels  at  RAE  Farnborough  for  a  number  of  years. 

2  DETAILS  OF  THE  DATA  INPUT  SCHEME 

The  scheme  provides  routines  for  reading  from  one  to  nine  real  variables, 
one-dimensional,  two-dimensional  and  three-dimensional  real  arrays,  and  similar 
routines  for  integer  variables  and  arrays.  The  scheme  is  based  on  a  routine 
FFORM  which  reads  in  a  line  of  data  as  individual  characters  and  then  builds  up 
the  appropriate  numerical  values  according  to  the  rules  proposed  in  section  2.1. 
In  sections  2.2  to  2.6  various  enhancements  to  this  basic  scheme  are  described. 


2 . 1  Number  recognition 

The  implementation  of  any  system  for  recognising  numbers  requires  a  defini¬ 
tion  of  a  number.  Such  a  definition  should  correspond  with  generally  accepted 
interpretations  and  be  as  simple  as  possible.  Here  it  is  defined  as  any  con¬ 
secutive  sequence  of  digits,  which  may  be  separated  by  a  limited  set  of  charac¬ 
ters  (digit  separators)  and  may  also  be  preceded  by  certain  other  characters 
(digit  prefixes).  The  digit  separators  are  a  decimal  point,  or  a  letter  E 
followed  immediately  by  a  digit,  a  space,  a  plus  sign  or  a  minus  sign.  The 
decimal  point  may  come  before  the  first  digit  (that  is,  be  part  of  the  prefix), 
and  the  characters  after  the  E  are  interpreted  as  an  exponent.  Any  number  may 
have  a  plus  or  minus  sign  as  its  first  character.  These  requirements  are  listed 

A  ... 

more  exactly  below  by  denoting  any  character  by  enclosing  it  m  brackets,  thus 
(>  .  For  example  the  character  <)  is  the  null  character  (that  is,  no  character), 
whereas  (  )  would  indicate  a  space.  Then  with  the  following  definitions, 


digit 

d 


E 

S 


(0>  , 

<1>  , 

(2)  ,  (3)  ,  <4>  , 

<5>  , 

(6)  ,  (7)  ,  (8)  , 

(9) 

(digit)  , 

d  (digit) 

( ie 

d  is  a  compact 

string  of 

digits) 

<•> 

<E>  , 

(E  > 

,  (E+>  ,  (E-) 

<  >  , 

<+>  , 

<-> 

the  numbers  must  be  in  one  of  the  eight  forms: 
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1  .  Sd 

2.  S.d 

3.  Sd.d 

4.  Sd. 

5.  SdEd 

6.  S.dEd 

7.  Sd .Ed 

8.  Sd.dEd 


In  the  Fortran  IV  routines  presented  here,  characters  are  read  sequentially  and 
a  number  is  detected  by  the  first  digit  found,  and  ended  when  the  character  list 
ceases  to  comply  with  one  of  the  forms  listed  above. 

2 . 2  Null  data 

A  common  requirement  in  inputting  data  is  that  an  input  variable  should 
retain  its  previous  value  (that  is,  a  null  datum  should  be  supplied).  In  the 
read  routines  described  here,  the  occurrence  of  two  commas  optionally  enclosing 
blanks  is  used  to  denote  a  null  datum.  Hence  the  character  line:  123,,  ,345 
will  give  four  values:  123,  two  null  data  and  345. 

A  comma  occurring  as  the  first  non-blank  character  on  a  line  implies  a  null 
datum,  but  a  comma  occurring  as  the  last  non-blank  character  on  a  line  does  not 
imply  an  extra  null  datum.  Hence  ,  123,345,  gives  three  values:  1  null 
datum,  123  and  345. 

It  should  be  noted  that  the  new  Fortran  77  free-format  specification 
differs  from  the  one  used  here,  in  that  Fortran  77  interprets  a  comma  occurring 
as  the  first  non-blank  character  on  a  line  as  a  null  datum  only  on  the  first 
line  of  data  input  by  a  read  statement.  Hence  two  lines  of  data: 

,123,345 

,567 

would  be  interpreted  by  the  input  routines  described  here  as  five  values:  1  null 
datum,  123,  345,  1  null  datum  and  567,  whereas  Fortran  77  would  interpret  these 
characters  as  four  values:  1  null  datum,  123,  345  and  567. 

2 . 3  Multiple  data 

Another  common  occurrence  in  inputting  data  is  that  several  consecutive 
items  in  the  data  are  identical.  In  this  case,  the  data  input  can  be  simplified 
by  using  the  convention  i*V  to  denote  i  consecutive  occurrences  of  the  value 
V  .  Specifically,  i  must  be  a  non-negative  integer.  Hence  3*3.14159  is 
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interpreted  as  three  values  of  3.14159  and  3*,  is  interpreted  as  three  null  data. 
When  used  in  this  way,  the  character  *  will  not  generate  a  warning  (see 
section  2.6).  Except  in  the  case  of  multiple  null  data  it  is  essential  that  the 
characters  of  i*V  should  be  compact.  For  example,  3*  3.14159  will  be 
interpreted  as  two  values  3  and  3.14159.  A  similar  interpretation  will  be  made 
of  3.0*3.14159,  or  3**3. 14159  or  3  *3.14159. 

If  this  convention  is  used  to  generate  more  data  than  is  required,  the 
input  routine  will  end  when  the  required  amount  of  data  has  been  read.  For 
example,  100*1.5  can  be  used  to  fill  a  60-element  array  with  data.  No  warning  is 
issued. 

2 . 4  Comments 

It  will  be  appreciated  that  the  definition  of  a  number  proposed  in  section 
2.1  permits  comments  not  containing  digits  to  be  included  in  the  data.  Hence 
for  example  123, ABC, 456  will  be  interpreted  as  123  and  456. 

As  the  facility  to  include  alphanumeric  comments  in  the  data  would  seem  to 
be  of  value,  a  convention  is  employed  in  which  characters  between  inverted 
commas  (either  single  or  double)  will  be  interpreted  as  a  comment.  Hence  the 
characters  _ _ „ 


'  AE3" 


will  be  interpreted  as  comments  but 


will  be  interpreted  as  the  value  6. 

More  specifically  an  alphanumeric  comment  is  preceded  by  a  <")  or  <’>  and 
ended  by  a  second  <">  or  <’)  or  the  end  of  a  line. 

2 . 5  End  of  read  characters 

While  the  user  must  specify  the  number  of  variables  to  be  read  by  the 
routines  (and  normally  the  routines  will  continue  reading  until  this  number  has 
been  found)  the  character  </)  will  end  the  read  call  forthwith,  any  further 
characters  beyond  </)  on  the  same  line,  being  ignored.  Thus  (/)  can  be  used  for 
the  protection  of  further  data  or  as  a  method  of  retaining  the  existing  values 
for  all  further  variables  of  the  read  call. 


TM  Ae  1804 


TM  Ae  1804 


7 


m 


For  all  the  routines  except  the  two-dimensional  and  three-dimensional 
array  routines,  $  has  the  same  effect  as  /.  For  two-dimensional  arrays,  $$  is 
needed  to  end  the  read  while  $  terminates  the  reading  for  the  current  row  or 
column.  Similarly  with  three-dimensional  arrays,  $$$  ends  the  read,  while  $  ends 
the  current  row  or  column  and  $$  ends  the  current  row  and  column. 

It  should  be  noted  that  $  and  /  have  no  effect  if  contained  within  an 
alphanumeric  comment  (see  section  2.4). 

2.6  Warnings 

Under  certain  conditions  the  routines  generate  warning  messages.  These 
messages  are  output  to  device  50,  but  can  be  suppressed  altogether  by  setting  the 
input  parameter  IDENT  to  zero.  It  should  be  noted  that  generation  of  a  warning 
will  not  stop  the  program.  However,  even  if  the  warning  messages  are  suppressed, 
the  output  parameter  IERR,  which  is  normally  positive,  is  set  negative  if  a 
warning  has  been  generated,  so  that  the  user  t..ay  take  appropriate  action  to  stop 
the  program  automatically  if  he  wishes. 

Details  of  the  warning  messages  are  given  in  section  3.2,  but  the  condi¬ 
tions  for  which  they  are  generated  are  listed  below.  In  the  input  scheme  des¬ 
cribed  here,  the  standard  data  delimiters  are  assumed  to  be  O  or  (,}  .  Other 
characters  can  appear  immediately  adjacent  to  data,  but  will  generate  warning 
messages  unless  otherwise  stated. 

Specifically,  if  a  number  does  not  start  at  the  beginning  of  a  new  line, 
is  not  immediately  preceded  by  (,),  <  )  ,  <">,  or  (')  or  is  not  immediately 
followed  by  <,),  <  ),  <->,  (+),  <"),  ('),  <$>,  </>  or  the  end  of  a  line,  a 

warning  message  locating  the  position  of  the  offending  character  is  generated. 

An  exception  is  the  character  (*)  when  used  to  denote  multiple  data  (see 
section  2.3).  Similarly  if  an  exponent  contains  other  than  (  >,  (+>,  (-)  or 
digits,  a  warning  is  generated. 

In  addition,  if  real  data  is  encountered  in  a  routine  for  reading  integer 
values,  the  value  is  set  to  the  nearest  whole  number  and  a  warning  is  generated. 

Examples  of  various  acceptable  data  lines  and  ones  which  will  generate 
warnings  are  given  below. 

(a)  123,456,789  or  123  456  789  will  r.ot  generate  a  warning; 

(b)  !23"ABC"456  will  not  generate  a  warning,  but  123ABC456  will  generate 
two  warnings; 


i 


m 
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(c)  +3.78-2.63-1.97+2.96  will  not  generate  a  warning; 

(d)  52*1.97  will  not  generate  a  warning  since  multiple  data  will  be 
assumed,  but  52**1.97  will  generate  two  warnings; 

(e)  1,3, A. 2, 5  ,  if  read  by  an  integer  read  routine  will  generate  a 
warning  and  the  third  value  will  be  set  to  A. 


3  THE  PROGRAMS 


The  programs  comprise  2A  subroutines: 


READ  1 ,READ2 , . . . , READ9 

READI 1 , READI2 , . . . .READI9 

READA 

READIA 

READB 

READIB 

RE  ADC 

READIC 


-  read  from  one  to  nine  real  variables, 

-  read  from  one  to  nine  integer  variables, 

-  read  a  one-dimensional  real  array, 

-  read  a  one-dimensional  integer  array, 

-  read  a  two-dimensional  real  array, 

-  read  a  two-dimensional  integer  array, 

-  read  a  three-dimensional  real  array, 

-  read  a  three-dimensional  integer  array. 


In  addition  the  following  subroutines  are  used: 


VAR9 

IVAR9 

ARRAY 


FFORM 


-  called  by  READI ,... ,READ9 .  VAR9  reads  nine 
real  variables  using  READA, 

-  called  by  READI 1 ,... ,READI9 .  IVAR9  reads  nine 
integer  variables  using  READIA, 

-  called  by  READA,  READB,  RE ADC,  READIA,  READIB, 
READIC.  ARRAY  reads  real  or  integer  variables 
into  an  array, 

-  called  by  ARRAY.  FFORM  translates  a  line  of 
characters  into  numerical  values  according  to 
the  rules  proposed  in  section  2. 


All  the  routines  are  written  in  standard  Fortran  IV.  The  characters  <"> 
and  <’>,  while  not  included  in  the  standard  character  set  are  available  in 
Fortran  on  most  computers  and  have  been  used  in  the  routines.  Within  the 
routines  it  is  assumed  that  characters  are  input  from  unit  IUNIT  and  that  warn¬ 
ings  are  printed  on  unit  50.  Each  "READ"  statement  should  be  associated  with  a 
whole  number  of  lines  of  data  (that  is,  data  for  two  "READ"  statements  should 
not  occur  on  the  same  line). 


In  the  remainder  of  this  section  the  routines  are  specified  in  more  detail. 
The  "READ"  routines,  together  with  VAR9  and  IVAR9  are  fairly  straightforward  in 
terms  of  logic  and  hence  these  are  presented  with  a  brief  specification  and 
listing.  Subroutines  ARRAY  and  FFORM  are  logically  more  complex  and  flow 
diagrams  and  annotated  listings  are  included. 
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3.1  Specification  of  subroutines 


3.1.1  READ1 ,READ2, . . . , READ9,READI1 >READI2 . READI9 


Purpose 


To  read  from  one  to  nine  real  (READn)  or  integer  (READln) 
variables  in  free  format. 


Call  statements 


CALL  READn (IUNIT, IDENT, IERR,X 1 , . . .Xn) 
or 

CALL  READln (IUNIT , IDENT ,IERR,L 1 , . . .Ln) 
where  1  <  n  <  9  . 


Input  parameters 


Output  parameters 


IUNIT  =  unit  identification  fcr  data  input 

IDENT  =  statement  identifier  and  warning  control: 

If  IDENT  =0,  no  warnings  are  printed; 

If  IDENT  S'  1 ,  identified  warnings  are  printed; 
If  IDENT  S'  1000,  all  data  are  printed  too. 

IERR  =  number  of  values  read  x  sgn 

where  sgn  =  1  if  there  are  no  warnings, 

=  -1  if  there  are  warnings. 

XI, ...Xn  or  Ll,...Ln  (1  ^  n  <  9)  are  set  to  the  values 
read  in. 


Subroutines  called  READn  uses  VAR9,  READA,  ARRAY  and  FFORM, 

READln  uses  IVAR9,  READIA,  ARRAY  and  FFORM. 

Exit  The  routines  end  when  either  n  values  have  been  read  in 

or  a  </)  or  a  ($>  is  encountered. 


3.1.2  READA,  READIA 


Purpose 


Call  statements 


To  read  a  one -dimensional  real  (READA)  or  integer  (READIA) 
array. 

CALL  RE ADA (IUNIT, IDENT, IERR,A, 1 1 ,12) 


Input  parameters 


CALL  READIA(IUNIT, IDENT, IERR, IA, 11 ,12) 

IUNIT  and  IDENT  are  defined  in  section  3.1.1; 

A  is  a  one-dimensional  real  array  of  dimension  S*I2, 

IA  is  a  one-dimensional  integer  array  of  dimension  ^12, 

II,  12  are  the  initial  and  final  subscripts  for  which  data 
is  read. 


Output  parameters 


IERR  is  defined  in  section  3.1.1. 

The  values  read  in  are  stored  in  A  or  IA. 


Subroutines  called  Both  routines  call  ARRAY  and  FFORM. 


The  routines  end  when  either  12— I 1+1  values  have  been  read 
in  or  a  (/)  or  a  ($>  is  encountered. 
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3.1.3  READB, 
Purpose 

Call  statements 

Input  parameters 


Output  parameters 

Subroutines  called 
Exi  t 


3.1. A  RE ADC, 
Purpose 

Call  statements 

Input  parameters 

\ 


READ IB 


To  read  a  two-dimensional  real  (READB)  or  integer  (READIB) 
array. 

CALL  READB (IUN IT, IDENT , IERR, B, 1 1 ,I2,NI,J1 ,J2,NJ,I0RD) 
or 

CALL  READIB (IUNIT , IDENT, IERR, IB, 1 1 ,12 ,NI , J 1 ,J2,NJ,I0RD) 

IUNIT  and  IDENT  are  defined  in  section  3.1.1, 

II,  12  are  the  initial  and  final  values  of  I,  and 
J 1 ,  J2  are  the  initial  and  final  values  of  J  for  which  data 
is  required  in  the  two-dimensional  arrays  B(I,J)  or  IB(I,J). 
NI,  NJ  are  the  dimensions  I  and  J  respectively  of  B  or  IB. 
IORD  denotes  the  order  in  which  the  array  is  filled.  If 
IORD  =  21,  the  subscript  I  varies  while  J  remains  fixed 
(columns  are  read  in  matrix  terms).  Otherwise  J  varies 
while  I  remains  fixed  (rows  are  read  in  matrix  terms) . 

IERR  is  defined  in  section  3.1.1. 

The  values  read  in  are  stored  in  B  or  IB. 

Both  routines  call  ARRAY  and  FFOEM. 

The  routines  end  when  (12- 1+ 1 )* (J2-J 1  + 1 )  values  have  been 
read  in,  or  when  a  (/>  or  two  consecutive  ($>  characters 
are  encountered.  In  these  routines,  a  single  ($>  denotes 
the  end  of  a  sequence  of  the  more  rapidly  varying  sub¬ 
script.  It  should  be  noted,  however,  that  the  appearance 
of  a  (/)  or  <$>  on  a  line  causes  further  data  on  that  line 
to  be  ignored.  Hence  two  sequences  terminated  by  ($) 
should  not  appear  on  the  same  line. 


READIC 


To  read  a  three-dimension  .  real  (READC)  or  integer 
(READIC)  array. 

CALL  READC (IUN IT, IDENT, IERR, C, II , 12 ,NI , J 1 , J2 ,NJ ,K1 ,K2,NK, 

I0RD) 

or 

CALL  READIC (IUNIT, IDENT, IERR, IC, II ,I2,NI,J1 ,J2,NJ,K1 ,K2,NK, 
I9RD) 

IUNIT  and  IDENT  are  defined  in  section  3.1.1, 

II,  12  are  the  initial  and  final  values  of  I, 

Jl,  J2  are  the  initial  and  final  values  of  J,  and 

K1 ,  K2  are  the  initial  and  final  values  of  K  for  which  data 

is  required  in  the  three-dimensional  arrays  C(I,J,K)  or 
IC(I,J,K).  NI,  NJ,  NK  are  the  dimensions  I,  J  and  K 
respectively  of  arrays  C  or  IC. 

I0RD  denotes  the  order  in  which  the  array  is  filled. 

I9RD  =  132  means  that  J  (identified  with  2)  varies  most 
rapidly,  followed  by  K  (identified  with  3),  followed  by  I 
(identified  with  1).  Other  significant  values  of  IORD 
are  as  follows: 
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I6RD  =  231  implies  JKI  order, 

=  213  implies  JIK  order, 

=  312  implies  KIJ  order, 

=  321  implies  KJI  order, 

For  all  other  values  of  I8RD,  the  order  is  assumed  to  be 

IJK. 

IERR  is  defined  in  section  3.1.1. 

The  values  read  in  are  stored  in  C  or  IC. 


Subroutines  called  Both  routines  call  ARRAY  and  FFORM. 


The  routines  end  when  (I2-I 1+1 )* (J2-J I +1 )* (K2-K1+ 1 )  values 
have  been  read  in  or  a  </>  or  three  consecutive  <$>  's 
are  encountered.  In  this  routine  a  single  ($}  denotes  the 
end  of  a  sequence  of  the  most  rapidly  varying  subscript  and 
two  consecutive  (?)  characters  denote  the  end  of  a  two- 
dimensional  sub-array  with  the  two  most  rapidly  varying 
subscripts . 


3.1.5  VAR9,  IVAR9  (subsidiary  routines) 


Purpose 


Call  statements 


Input  parameters 


To  read  N  real  (VAR9)  or  integer  (IVAR9)  variables. 

CALL  VAR9 (IUNIT,IDENT,IERR,N,X1 , . . . ,X9) 
or 

CALL  IVAR9 ( I UN IT  ,IDENT , IERR , N , L 1 , . . . , L9 ) 

IUNIT  and  IDENT  are  defined  in  section  3.1.1. 

N  is  the  number  of  variables  required  1  <  N  <  9  . 


Output  parameters  IERR  is  defined  in  section  3.1.1. 

XI,... XN  or  L1,...,LN  are  the  values  found 

Subroutines  called  VAR9  uses  READA,  ARRAY  and  FFORM, 

IVAR9  uses  READIA,  ARRAY  and  FFORM. 

Exit  The  routines  end  when  either  N  values  have  been  read  in  or 

a  </>  or  a  <$)•,  is  encountered. 


3.1.6  ARRAY  (subsidiary  routine) 


Purpose 


Call  statement 


Input  parameters 


To  read  data  into  a  three-dimensional  real  or  integer  array. 

CALL  ARRAY (IUNIT, IDENT, IERR, C,NI,NJ,NK,IC, MI, MJ,MK,LIM, 
J6RD,IDIM,ININT) 

IUNIT,  IDENT  are  defined  in  section  3.1.1 

ININT  =  1  if  ARRAY  is  called  by  an  integer  "READ"  routine 
=  0  otherwise 

NI,  NJ,  NK  are  the  dimensions  of  a  real  array  C(I,J,K) 

MI,  MJ,  MK  are  the  dimensions  of  an  integer  array  IC(I,J,K) 

IDIM  is  the  dimension  of  array  in  "READ"  statement 
(=1  for  reading  variables) 
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Output  parameters 

Subroutines  called 
Exit 

An  annotated  listing 
in  Appendix  A. 

3.1.7  FFORM  ( 

Purpose 

Call  statement 
Input  parameters 


Output  parameters 


Subroutines  called 
Exit 

An  annotated  listing 
in  Appendix  A. 


LIM(6)  is  an  array  containing  initial  and  final  values  of 
I,  J,  K  for  which  data  is  required 

J9RD  indicates  the  order  in  which  the  array  is  filled  (see 
section  3. I . A) . 

IERR  is  defined  in  section  3.1.1. 

The  values  read  in  are  stored  in  IC  if  ININT  =  1  or  C 
otherwise . 

FFORM 

The  routine  ends  when  the  required  number  of  data  has  been 
read,  when  (/>  is  tound  or  when  ($>  is  encountered  for 
IDIM  =  1,  ($$>  is  found  for  IDIM  =  2,  ($$$>  is  found  for 

IDIM  =3. 

of  ARRAY  appears  in  Appendix  B  and  a  flow  diagram  is  given 


subsidiary  routine) 


To  read  a  line  of  characters  and  translate  into  numerical 
values . 

CALL  FFORM (IUN IT, IDENT ,NREC, INEND, INWAR,NUM,AA, JDAT, ININT) 

IUNIT  and  IDENT  are  defined  in  section  3.1.1 
NREC  is  the  data  record  counter 

INWAR  is  the  number  of  warnings  generated  in  current  read 
routine 

ININT  =  1  for  reading  integers 
=  0  otherwise 

AA(40)  is  an  array  for  storage  of  numbers  found 

INEND  =  0  if  no  end  of  read  character  has  been  found 

1  if  ($>  has  occurred 

2  if  <$$)  has  occurred 

3  if  <$$$>  has  occurred 

4  if  (/)  has  occurred. 

INWAR  is  the  number  of  warnings  generated  in  current  read 
routine 

JDAT(4o)  is  an  array  to  indicate  if  the  corresponding  value 
in  AA  is  normal  (JDAT  =  1),  null  (=  0)  or  multiple  (=  2) 

NUM  is  the  number  of  values  found. 

None 

The  routine  ends  when  the  end  of  a  line  is  reached, 
of  FFORM  is  given  in  Appendix  B  and  a  flow  diagram  appears 
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3.2  Warning  messages 

The  conditions  for  which  warning  messages  are  generated  are  defined  in 
section  2.6.  The  warning  messages  are  normally  printed  on  device  50,  but  can  be 
suppressed  by  setting  the  parameter  IDENT  to  zero. 

The  first  warning  condition  found  in  a  routine  gives  the  following  heading: 

****  DATA  READ  WARNINGS  **** 

Then  one  or  more  of  the  following  warning  messages  will  be  given: 

(a)  FOR  IDENT  =  I,  AFTER  J  LINES  &  K  CHAR  C^  message, 

where  message  is  either  STARTS  NUMBER  , 

or  ENDS  NUMBER  , 

or  IN  EXPONENT  . 

Here  I  =  identifier 

J  =  number  of  lines  of  data  read  in 
K  =  character  count  on  current  line 
CjC^  are  the  two  characters  in  positions  K-l  and  K  . 

(b)  FOR  IDENT  =  I  ,  V  FOUND,  INTEGER  N  ASSUMED, 
where  I  =  identifier 

V  =  floating  point  value 
N  =  nearest  integer  value  to  V  . 

If  IDENT  ^  1000  and  a  warning  message  has  been  given,  the  line  of  characters  is 

printed  out  followed  by  a  blank  line  with  the  character  +  showing  the  positions 

of  the  characters  generating  any  warnings. 

3.3  Implementation 

The  free-format  data  input  scheme  presented  here  has  been  implemented  and 
tested  on  a  number  of  computer  systems.  Up  until  now  the  routines  have  been 
used  successfully  on  the  following  computers: 


I) 

DEC  KL10B 

at 

RAE  Farnborough, 

2) 

IBM  370/168 

at 

UKAEA  Harwell, 

3) 

ICL  1906S 

at 

RAE  Farnborough, 

4) 

ICL  4130 

at 

RAE  Bedford, 

5) 

Honeywell  Mk 

III 

'  Time  Sharing  Service. 

The  only  modification  found  to  be  necessary  was  on  the  ICL  4130,  where  the  double 
quotation  character  <">  is  not  allowed  in  Fortran  programs.  For  use  on  the  4130, 
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therefore,  all  references  to  the  variable  IDQ  in  routine  FFORM  should  be  deleted. 
(The  inclusion  of  alphanumeric  comments  can  still  be  achieved,  of  course,  by  the 
use  of  single  quotes  (') .) 

The  core  requirements  for  the  data  input  scheme  depend  on  how  many  of  the 
'READ'  routines  are  called.  For  example,  each  of  the  individual  READn  routines 
requires  approximately  40  words,  but  the  call  will  also  involve  READA  (80  words), 
VAR9  (200  words),  ARRAY  (600  words)  and  FFORM  (1800  words).  Calling  all  24 
'READ'  routines  would  need  approximately  4K,  but,  for  typical  usage,  involving 
say  six  of  the  routines  the  storage  requirement  would  amount  to  3K. 

Experience  has  shown  that  the  'READ'  routines  are  between  10  and  15  times 
slower  than  using  the  normal  Fortran  formatted  input.  Typically,  on  the 
DEC  KUO,  1000  numbers  can  be  input  in  less  than  10  seconds. 

On  the  RAE  DEC  KL10B,  the  read  system  is  implemented  for  the  users’ 
convenience  in  SYS:  RAE  LIB.  A  help  file  HELP  F0READ  gives  details,  and  a 
program  called  by  RUN  READEM  gives  a  demonstration  of  the  system.  In  practical 
use  on  the  DEC,  the  read  system  has  been  found  to  be  reliable  and  very 
accommodating.  The  error  warnings  on  the  DEC  are  arranged  to  appear  on  device 
5,  that  is  at  the  user's  teletype  or  in  the  LOG  file  for  batch  runs. 

4  CONCLUDING  REMARKS 

Routines,  written  in  Fortran  IV,  for  the  input  of  numerical  data  in  free- 
format  have  been  described.  Provided  a  set  of  simple  rules  are  obeyed,  the 
scheme  gives  the  user  the  flexibility  to  format  data  in  the  most  convenient 
way.  Additional  features  include  provision  for  input  variables  to  retain  their 
previous  values  (null  data),  for  identical  consecutive  data  to  be  described  in 
a  simplified  form  (multiple  data)  and  for  alphanumeric  comments  to  appear 
amongst  the  data. 

The  input  scheme  was  devised  primarily  to  simplify  the  reading  of  data 
into  large  flow-field  calculation  programs  which  are  being  developed  at  RAE  for 
use  throughout  the  UK  aircraft  industry.  For  this  type  of  program,  the  extra 
core  requirements  and  time  taken  when  reading  in  data  in  this  way  (typically 
less  than  4K  and  less  than  10  seconds  per  1000  numbers  on  the  DEC  KL10B)  are 
very  small  compared  with  the  resources  used  by  the  program  as  a  whole  (typically 
•b4K  and  -30  minutes  running  time).  The  use  of  standard  Fortran  IV  throughout 
the  routines  means  that  programs  using  this  input  scheme  can  be  transferred 
easily  from  one  computer  to  another. 
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Appendix  A 
FLOW  DIAGRAMS 
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Appendix  A 


FLOW  DIAGRAM  OF  SUBROUTINE  FFORM  (2) 
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Appendix  B 
LISTINGS 


0031 
0032 
0333 
0034 
0035 
0036 
0033 
6333 
0039 
0:j  i  ii 
00  1 1 
00  ■  2 
6013 

00  15 
00 . 6 
0017 
00  IS 
0019 
0-j23 
0.J21 
0022 
f  3  23 
0324 
00.25 
0b26 
00  4  7 
0020 
03  2? 
04  23 
0U31 
C03.? 
0333 
0034 
0335 
0  b  25 

%u?y 

0>  ■  55 
0U59 

•jj  -.3 

W*  1 
&34J 
004?, 
0i  144 
0..  15 
00  45 


30-0 

ro ' ; 

0ij33' 
C :  .  4 


C>MotoK*>to(c4oK*>M<***  **OMoK*********>Mc:to|OK*:  **>MOK*)K*¥***:f. *******  **>tofc**3^****¥*** 


C 

C  FORTRAN  FRF.E-FCRMAT  DATA  INPUT  SCHEME 

C 

C  1  :  RE ADI  TO  READS  &  READ  1 1  TO  READ!  9 

C 

C  G  SUTLER  &  J  PIKE  :  MAr'  T3 

C 

C*#***#****»********#**:<>:*¥****W**  ********  i**********^**************** 

c 

C  THESE  ROUT INFS  ALLOW  FROM  1  TO  S  PEAL  OR  INTEGER 

C  '7APIAE.EE  TO  RE  REaD 

C 

c 

C  INPUT  PARAMETERS 

C  - 

C 

C  I DENT  =  IDENTIFIER  AMD  WARNING  CONTROL 

C  IF  I DENT  -  g.  NO  WARNINGS  ARE  PRINTED 

L  IF  IDtiNT.Gc.  1,  IDENTIFIED  WARNINGS  ARE  PRINTED  ON 

L  DEVICE  50 

C  IF  iDtNT.GE. 1000,  ALL  DATA  ARE  PRINTED  ON  DEVICE  50 

L  I  UN  IT  =■-  DATA  i.iFUT  UNIT 

C 

c 

C  OUTPUT  PARAMETERS 

C  - 

c 

0  !EPPs  =  ERROR  INDICATOR 

C  =  NUM.-CR  0 f  VALUES  READ  MULTiFL JED  Bv  SGN 

C  tJHE.'c  SGli  «  1  IF  NO  WARNINGS  HAVE  OCCURREAD 

C  SGN  =  -1  IF  WARNINGS  HAVE  OCCURkEAi) 

C  XI. . XM  OR  LI, . L!1  ( 1  .LE.M.LE.Sl  ARE  SET  TO  THE 

C  VALUES  FOUND 

C 

C 


C»>i<>|olo,«))o)o)(*)lcK*¥:::¥i.':;>-¥*¥iK>i<:i<.')(¥¥*4:¥¥>i¥¥.¥4'^,k.i<¥¥:S¥*‘«¥:.~<;!<\!*:.'R)|.  b1;*  **************** 


L 


SUBROUTINE  READ  1  ■.  JUnIT,  IDENT.  I  ERR, XI) 

CALL  VAR9UUNIT.  IDENT,  I  ERR,  l,X.,,;2,X3»X4,X5,XS.X?,X8,X9) 

RETURN 

END 

SUE ROUTINE  READ? (I UN  IT, IDENT, IERR.X1.X2/ 

CALL  VAR  3  I  UN  iT,  'TENT,  I E  k  k ,  2 ,  X ". ,  X2 , ;  3 ,  :{A .  X5 ,  XS ,  XT ,  X? .  X9 ) 

Rs’uRH 

END 

SU'.’i'.OUTINE  READS  ■! I  CM  3  T,  1 1  ENT,  I  ERR ,  X 1 ,  X2 ,  X3 ) 

Li  13.  ’/ARP t  IUI1  IT,  IDENT,  iERr  ,3,’  II  ,X2.>G, X4.X5.X3. X?,  X3.X9) 
PHI.  AN 


::  NX 

CLLLCUTI  p*r.M(  I  UNIT, 
1  r  ...  T  :■  1 1  HIT,  TLENT,  I 

i.  .  i 


IDENT, ! ERR, Xi.X2.X3.  <0; 

-  Hk  .4.  XI.  X2.X3 ,  X4,  ,  X  J ,  X? , 


XP.X9) 
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0055 
0056 
005? 
0050 
0059 
0O50 
0051 
0052 
0653 
0054 
0055 
0056 
0057 
0058 
0059 
0073 
0071 
0072 
0073 
0074 
0075 
0076  C 
0077 
0073 
0079 
0030 
0031 
0032 
0033 
6054 
0635 
0636 
0037 
0038 
0039 
0090 
0031 
0092 
0093 
0094 
0035 
0036 
0697 
0098 
0099 
0100 
0101 
0132 
0103 
0134 
0135 
0106 
0107 
0108 


El'*-- 

SUBROUTINE  READS  ( IUN  IT,  IDENT,  IERR,  XI ,  X2, X3,  X4,  X5) 

CALL  VARS (IUNIT. IDENT. IERR, 5 .  X! .  X2.  X3.  X4.  X5.  X6,  X? .  X8.  X9) 

RETURN 

END 

SUBROUTINE  REAB6C I  UNIT,  IDENT.  IEPR.X1.X2,X3.X4,X5,XS) 

CALL  V4R9'  IUN  IT..  IDENT.  !ERP,S,Xi,X2,X3,X4,X5,Xb.X7,X8,X9) 

RETURN 

END 

SUBROUTINE  READ?( IUN IT, IDENT, IERR,X1,X2,X3.X4,X5.X5.X7) 

CALL  VAR9  ( IUN  IT,  I  DENT,  1  ERR .  7 .  X  >. ,  X2 ,  X3 ,  X4,  X5 .  X6 .  X7 ,  X8 ,  X9) 

RETURN 

END 

SUBROUTINE  READ8( IUNIT, IDENT, !ERR.X1.X2.X3,X4.><5.XS,X?.X8> 
COLL  VARS (IUN IT, IDENT, IERR.8,Xi»X2,X3.  X4,X5,X6,>;7,X8,X9) 
RETURN 
END 

SUBROUTINE  READ9 (IUNIT, IDENT, I ERR , X 1 , X2 , X3 , X4 . X5 , X5 . X7 , X8 , X9 ) 
CALL  VAR!?(  IUNIT,  IDENT.  IERR,9,X:,X2,X3,  X4,X5,X6,X7,X8,X9) 
RETURN 
END 

SUBROUTINE  READI 1 ( IUNIT, IDENT. IERR,L1) 

CALL  I VARS ( IUNIT, IDENT, IERR. 1,L1.L2,L3,L4,L5.LS,L7,LS,L9) 

RETURN 

END 

SUBROUTINE  READI2C  I'JNIT,  IDENT,  IERR,L1,L2) 

CALL  IVAR9( IUNIT, IDENT, IERR, 2, L 1 , L2, Lo, L4,L5,Lo,L7,L8,L9) 

RETURN 

END 

SUBROUTINE  READI3C IUNIT, IDENT, 1ERR.L 1 ,L2,L?) 

CALL  IVAR9C IUNIT  IDENT, IERR, 3. _ 1,L2, L3, L4,L5, L6. L7. L8,L9) 

RET URN 
END 

SUBROUTINE  READ I4( IUN IT, IDENT, IERR.Ll,L2.L3.L4> 

ChlL  I  VAR  3  (IJJNIT,  IDENT,  IERR, 4,  J  ,L2,L?  .1.4.  L.5  . 1 . 17,  L8.L9) 

B  E  (  URN 
END 

SUBROUTINE  READ  15 ( IUNIT, IDENT. IERR,L1.L2,L3,L4,L5) 

CALL  I  VAR  9  (  IUN  IT,  IDENT.  IERR, 5,  L 1 ,  L2,  LI!,  1.4.  L5,  LG .  L7,  L8,  L9) 

RETURN 

END 

SUBROUTINE  READI6C IUNIT, IDENT, IERR.Ll.L2.L3.L4.L5.L6) 

CALL  IVORS!  IUNIT,  IDENT,  IEC'K.6,  •_  1  .L2,  L3,  L4.L5,  Lb  .  L7,  LS,  L9) 

RETURN 

END 

SUBROUTINE  READ  17 ( IUNIT, IDENT, IERR,L1,L2,L3,L4.L5._6,L?) 

CALL  IVAR9( IUNIT, IDENT, IERR,7,L1,L2,L3.L4,LS,LS,L7.L8,L9) 

RETURN 

END 

SUBROUTINE  READ  13 C IUNIT, IDENT. IERR, LI .L2,l3,L4,L5,L6,L7,L8) 
CALL  IVAR9C IUNIT. IDENT, IERR.8,L1,L2,L3,L4,L5,L6,L7,L8,L9) 
RETURN 
END 


o 

00 


01 

< 


0139  SUBROUTINE  READ  19 ( IUN IT, IDENT, IERR,L1,L2,L3.L4.L5,L6,L7,L8,L9) 

0110  CALL  IVAR9'  IUNIT,  IDENT,  IERr,,9,_i,L2,L3.L4,L5,L6,L7.L8,L9) 

0111  RETURN 

0112  i-NI' 

****  LIST  END 


x 

H 
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0lj  3  X 

0032  C 

0333  C  FORTRAN  FREE-FCRMAT  BATA  INPUT  SCHEME 

0034  C 

0O35  C  2  :  READA,  READ  I  A.,  READB,  PEADIB,  READC  fc  READIC 

333b  C 

003?  C  G  BUTLER  8,  J  PIKE  :  MAY  '73 

0030  C 

0310  C 

3011  C  REfiDA  READS  A  ONE-DIMENSIONAL  REAL  ARRAY 

0012  C  READB  READS  A  TUO-D IMENS IONAL  REAL  ARRAY 

0013  „  READC  READS  A  THREE-DIMENSIONAL  REAL  ARRAY 

6314  C  READIA  READS  A  ONE-DIMENSIONAL  INTEGER  ARRAY 

0015  C  READ  IB  READS  A  TUG-D IMENS 1 3'<AL  INTEGER  ARRAY 

03  IS  C  READ  IE  READS  A  THREE-DIMENSIONAL  INTEGER  ARRAY 

30  1 7  C 
0913  C 

0919  C  INPUT  PARAMETERS 

P  320  C  - 

032 1  C 

0922  C  I DENT  =  IDENTIFIER  AND  WARNING  CONTROL 

0023  C  I UN  IT  =  DATA  INPUT  UNIT 

0024  C  A  IS  P  GNF.-P IMENS  1CNAL  REAL  ARRAY  OF  DIMENSION  .GE.I2 

001"  C  0  13  A  TWO-DIMENSIONAL  REA.  ARRAY  OF  DIMENSION  NI,NJ 

0G26  0  C  IS  A  THREE-DIMENSIONAL  REAL  ARRAY  OF  DIMENSION  NI,NJ,NK 

6027  C  IA  IS  «  OHE-P IMENS IONFL  INTEGER  ARRAY  OF  DIMENSION  .GE.I2 

0023  C  10  IS  A  TWO-DIMENSIONAL  INTEGER  ARRAY  OF  DIMENSION  NI.NJ 

0929  C  IC  IS  A  THPEE-P IMENS IONAL  INTEGER  ARRAY  OF  DIMENSION  NI,NJ,NK 

0020  C  11.12  ARE  THE  INITIAL  AND  -INAL  VALUES  OF  !, 

0031  C  J1.J2  ARE  THE  INITIAL  AND  -INAL  VALUES  OF  J, 

0032  ■:  K1..K2  ARE  THE  INITIAL  AND  -INAL  VALUES  OF  K, 

euaS  C  'OR  WHICH  DATA  IS  TO  3E  READ 

0534  C  IORD  DENOTES  THE  ORDER  IN  WHICH  THE  APPROPRIATE  ARRAY  IS  FILLED 

0235  C 
0036  C 

0137  C  OUTPUT  PARAMETERS 

6033  C  - 

iiA  35  C 

00 -'O  C  1ERR  =  ERROR  INDICATOR 

6041  C  TEE  VALUES  ARE  STORED  IN  THE  APPROPRIATE  ARRAY 

EC-T2  C 
0043  c 

0045  C 
Fa  43  C 

004?  SUBROUTINE  READAC I'.M'IT, IDENT, I  ERR, A, 1 1, 12 J 

0043  DIMENSION  AC  12  • ,  I  ft  *.  1)  ,LiM<.6) 

6049  DO  1  K  =  1 ,  G 

0053  1  Li'MC!0»l 

03-:  1  L  I  i'll  3 )  =  1 1 

L.352  L  111(6)  - 12 

00 53  JOPD- !?3 

3054  CALL  ARRAYCIUNIT, IDENT. 1EPR.A, 1.1, 12, IA, 1. 1, l.LIM.JORD. 1,0) 
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0055  RETURN 

0056  END 

005?  C 

0058  SUBROUTINE  READ I ft (I UN IT, IDENT, IERR, IA, II, 12) 

0059  DIMENSION  IAH2) ,  AC  1)  ,LIM(6) 

0050  DO  i  K=l,6 

0051  1  L IMCK.)  = 1 

0052  L IMC3)  =  1 1 

0053  LIMcS) =12 

0064  JORD= 123 

6055  CALL  ARRAYC IUNIT,  IDENT,  IERR, A, 1,1,1, IA, !, 1. 12.LIM, JORD, 1,1) 

0856  RETURN 

0057  END 

0053  C 

0059  SUBROUTINE  READS'!  IUH  IT,  IDENT,  IERR,  B,  II,  12, NI,  Jl.  J2,NJ,  IORD) 

0070  DIMENSION  BCNI.NJ) , !BC 1) ,L1MC6) 

0071  LIMC !)  =1 

0072  LIM(2)=I1 

0073  LIMCS) =J! 

0074  LIM<4) = l 

0075  L 1 M  <  5 )  =  1 2 

0076  LIMCS) =J2 

067?  JORD= 123 

0073  IFCiORP.EQ.21)  J0RD=132 

0079  CALL  ARRA 7 C IUNIT,  IDENT. IERR,B,  1,NI,NJ, IB, 1, 1, l.LIM, JCRD,2,0) 

0030  RETURN 

0031  END 

0082  C 

0033  SUBROUTINE  READ  IB  HUN  IT, IDENT, IERR, IB, 1 1 , 12.NI , J 1 , J2, NJ, IORD) 

0634  DIMENSION  I3CN1.NJ) ,BC 1) ,LIMc6) 

0035  LIMC l) =1 

0036  LIMC2)  =  1 1 

003?  L IMC?) =J 1 

0038  L  if1C4>  =  1 

0039  L  IM',5)  =12 

0093  L1IKS)  =J2 

0U91  JORD= 123 

0092  IFCiORD.EQ.21)  JCRD=132 

0093  CALL  RRRAYC  JUNIT, IDENT, IERR, B, 1, 1, 1, IB, l.NI,NJ,LIM. JORD, 2, 1) 

0094  RETURN 

0095  END 

0096  C 

009?  SUBROUTINE  READCC il'NIT, IDENT, IERR.C, I !, 12, NI. Jl, J2.NJ, 

005S  1  K1,K2,NK, IORD) 

0099  DIMENSION  CCNI.NJ.NK) , ICC  1) ,LI4C6) 

0130  LIMC I ) =11 

0131  L IMC2) =J 1 

0132  L IMC3) -K 1 

0103  L IMC4) = 12 

0134  LIMC5;=J2 

0105  L  IMCf..’  =K2 

0196  JCr.D=  12: 

0137  IrC  IORD. En.  132)  J0C'D=132 

0153  lFCinfD.ECi.23n  JORD =231 
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3;33  If-'(  I0RD.EQ.213)  J0RD=213 

0110  IFiI0RD.FQ.312)  J0RD=2 12 

0111  IF( I0RD.EQ.321)  J0RP-321 

ill?  CALL  ARRAYC IUNIT,  IDENT,  lERR.C.NI.NJ.NK,  IC,  1, 1.  l.LM,  JORD.3,0) 

3113  RETURN 

01 14  END 

01.5  C 

PUS  SUBROUTINE  READ  ICC  IUNIT,  IDENT,  I  ERR,  IC, II. I2.NI, Jl, J2.NJ, 

011?  1  K1,K2,NK, IORD) 

0118  DIMENSION  IC(NI,NJ,NK) ,LIM(6) ,C( 1) 

0119  LIM( 11=11 

0120  L 1M(2) =J 1 

&121  LIW3}=K’. 

0122  L  Iil(-i)  =  12 

0123  L IMC5J  =J2 

6124  LIM(6l=K2 

0125  J0RD-123 

0126  IFnORD.EQ.  132)  J0RD=132 

012?  IF ( I0RD.EQ.231)  JORD-231 

0128  IF ( I0RD.EQ.213)  JORD-213 

0129  IF( T0RD.EQ.312)  JORD-312 

0129  IF ( I0RD.EQ.321 )  JORD-321 

0131  CALL  HRRAYCIUNIT,  IDENT, IERR.C, 1,1,1. IC.NI.NJ.NK.LIM, J0RD.3, 1) 

0132  RETURN 

0133  END 

jkk#*  LIST  END  **** 
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000 1  C^W>Wor^)K>K>MoK*1oM<***^MoMo(oMoMo»M'**<!MoK*^*5MoMotoMoK*>toK*^*H<>loto(oWo(oto)oloK*>|o({ 

0002  C 

0003  C  FORTRAN  FREE-FORMAT  DATA  INPUT  SCHEME 

0034  C 

0035  C  3  :  VAR 9  8,  IVAR9 

8036  C 

003?  C  G  BUTLER  C,  J  PIKE  :  MAY  '?3 

0838  C 

0810  c 

0011  C  VARS  READS  N  REAL  VARIABLES  ( 1 .LE.N.LE.9) 

0012  C  IVAR9  READS  N  INTEGER  VARIABLES  Cl. LE.N.LE.9) 

ecus  c 

0314  C 

0015  C  INPUT  PARAMETERS 

00 i 6  C  - 

6017  C 

0013  C  IDENT  =  IDENTIFIER  a  WARNING  CONTROL 

0619  C  I UNIT  =  DATA  INPUT  UNIT 

0020  C  N  =  NUMBER  OF  VARIABLES  REQUIRED 

0021  C 

0022  C 

0023  C  OUTPUT  PARAMETERS 

0024  C  - 

0025  C 

032c  C  (ERR  =  ERROR-INDICATOR 

0027  C  XI, X2 . ,XN  OR  L1.L2 . LN  ARE  THE  VALUES  FOUND 

0028  C  UHERE  Cl. LE.N.LE.9) 

9023  C 

0030  C 

0031 
0032  C 

0033  C 

0034  SUBROUTINE  VAR9C I UN  IT, IDENT, IERR,N,X1,X2,X3.X4.X5.X6 

0035  1  ,X?,X3,X9) 

0636  DIMENSION  AH (9) 

0037  AA ( 1 ) =X1 

0038  AAC2) =X2 

0039  AA(3) =X3 

O040  ftH ( 4) - X4 

0041  AAC5) =X5 

0042  AA(6)=X6 

0043  AA  ■'  7 )  =X7 

0044  HA (8) =X8 

0045  AA  (9)  -=X3 

0046  CALL  REHDAUUNIT,  IDENT,  IERR.AP,  1,N) 

0847  XI-kAU) 

0048  X2=AA (2) 

0349  X3-AAC3) 

0050  X4-AA (4) 

0051  X5=AA(5) 

0052  X6  =Hh  i.  6 ) 

0053  X7=AAf"'' 

0054  X3=PAC8) 
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0055 

0056 

0057 

0053 

0059 

0050 

0051 

0062 

0053 

0054 

O055 

005b 

0057 

0058 

0053 

0070 

0071 

0072 

GG73 

0074 

0075 

0076 

0077 

0079 

0979 

0830 

0031 

0032 

XiXOjOfc 


1 


X3=AA<9) 

RETURN 

END 

C 

SUBROUTINE  IVAR9( IUNIT, IDENT. IERR,N,L1,L2,L3,L4,L5,L6 
1  ,L?.L3.L9) 

DIMENSION  LAC9) 

LA( 1) =L1 
LA  (2) =L2 
LA  (3) =L3 
LA  (4) =L4 
LA (5) =L5 
LfKS)  ~L6 
LA  (?) =L 7 
LA'.  9)  =L8 
LA (31 =L3 

CALL  RED1AC IUNIT. IDENT, IERR,LA, 1,N) 

L 1 =LA (1 ) 

l:?=la(2) 

L3=LA (3) 

L4=L0(4) 

L5=LA (5) 

L6=LA(6) 

L7  =LA (7) 

L3-LA(3) 

L9-LA  ('20 

RETURN 

END 

LIST  END  Xdotcfc 
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0031 

0032 

C 

0033 

C 

FORTRAN  FREE-FORMAT  DATA  INPUT  SCHEME 

0034 

C 

0035 

c 

4  :  ARRAY 

0036 

■  c 

003? 

c 

G  BUTLER  &  J  PIKE  :  MAY  '?3 

0033 

c 

0039 

0010 

C 

0011 

C 

ARRAY  READS  DATA  INTO  THE  THREE-DIMENSIONAL  REAL 

ARRAY  C  OR 

0012 

c 

THE  THREE-DIMENSIONAL  INTEGER  ARRAY  IC  ,  ACCORDING  TO  THE 

0013 

c 

VALUE  OF  THE  PARAMETER  IN I NT 

0014 

c 

0U15 

c 

0016 

c 

INPUT  PARAMETERS 

901? 

c 

— 

0013 

c 

00 13 

c 

IDENT  =  IDENTIFIER  &  URPNING  CONTROL 

0920 

c 

IUNIT  =  DATA  INPUT  UNIT 

01121 

c 

N I , NJ, NK  =  DIMENSIONS  OF  REAL  ARRAY  C 

0022 

c 

-  DIMENSIONS  OF  INTEGER  ARRAY  IC 

0023 

c 

I D I M  =  DIMENSION  OF  ARRAY  IN  CALLING  ROUTINE 

0324 

c 

ININT  =  1  FOR  READING  INTEGERS  ,  0  OTHERWISE 

0025 

c 

L I Ml 61  =  ARRAV  CONTAINING  INITIAL  AND  FINAL 

0326 

c 

VALUES  OF  !,J1<  FOR  WHICH  DATA  IS 

REQUIRED 

002? 

r 

w 

JORD  =  ORDER  iN  WHICH  the  ARRAY  IS  FILLED 

0023 

c 

0029 

c 

0030 

c 

OUTPUT  PARAMETERS 

0071 

c 

0032 

c 

0033 

c 

IERR=  ERROR  INDICATOR 

0034 

c 

VALUES  READ  IN  ARE  STORED  IN  IC  IF  iNINT=  1, 

0035 

c 

OR  C  OTHERWISE. 

0036 

c 

003? 

c 

G338 

c 

OTHER  PARAMETERS 

0039 

c 

0040 

c 

0041 

c 

AA (401  -  ARRAY  FOR  NUM3ER  RETURN  FROM  FFORM 

0342 

c 

INEND  =  ENDING  INDICATOR 

60  43 

c 

INLJAR  =  WARNING  INDICATOR 

0044 

c 

JDAT143)  =  ARRAY  TO  INDICATE  IF  CORRESPONDING 

AA 

VALUE 

EC  45 

c 

IS  HULL  OR  MULTIPLE 

00'S 

r 

KM  =  MULTIPLE  DAI  A  COUNTER 

004? 

c 

KN  =  COUNTER  FOR  STORING  CURRENT  LINE 

0043 

c 

MULT  =  MULTIPLIER  FOR  MULTIPLE  DATA 

0049 

c 

NCOUNT  =  COUNTER  FOR  NUMBERS  FOUND 

0038 

c 

NREC  =  DATA  LiH.1  COUNTER 

003! 

£ 

NUM  -  NUMBER  OF  VALUES  FPU'  D  ON  CURRENT  LINE 

3952 

c 

0053 

c 

6054 

r 

o 

00 
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Of  55 
0056 
6057 
0038 
0658 
0050 
805? 
0652 
0053 
0054 
0035 
0058 
06  5  7 
0058 
005.3 
6070 

C  ICOUNT,  JCOUNT,  KC0UN7  REFER  Tu  A  *  NfiTUffft'.  M 

C  COUNTING  ORDER,  UITH  K  VARYING  HOPE  RAPIDLY 

C  THAN  J,  WHICH  IN  TURN  VARIES  MORE  RAPIDLY  THAN 

C  I.  THE  I COUNT.  JCOUNT,  COUNT  V-»i_UES  "ji-  ANY 

C  DATA  ITEM  APE  THEN  RELATED  TC  ;V5  ADDRESS  IN 

C  ARRAY  C  OP  ?C  USING  II.  JJ,  KK  DERIVED  PROM  THE 

C  ORDER  PARAMETER  IORD. 

C 

L»-<J|o(olololo(o|o|oloh'.X:.:'i:Ho(cj(t**5K>i.''::!k«)KH<H:;*',<*H:>l;!K'K**>K)('*.:  ***  ****tar.)t  ****HoK#**>tc** 

C 

C 

SUBROUTINE  ARRAY ( IUNIT.  iDENT. IERR . C.  N I .  NJ, NK,  iC.MI.MJ.MK, 

1  LIM, JOED. IPIn. Im INTI 

DIMENSION  : f N I .. NJ, NK; ,  IC  CMI.MJ.HK)  ^0) ,  *LA7 r  43?  , LIM(6) 

C .  .  _ 

0071 

C 

INTERPRET  ORDER  AND  LIMIT  VALUES 

0672 

1 1  --J0RD/1Q0 

0073 

jj=t:joRD-;30*ii)/i0 

0074 

KK=JORD-!ed*n-ie*jj 

G675 

N 1 1 =L iM( i  1 J 

0076 

NI2=Lim  I  I+3) 

0077 

NJ  1  =L  IM(  J  I, 

0078 

IIJ2=L  IM(  JJ+3) 

807  J 

NK  1  --LIMCKIO 

0030 

NK2=L IMCKk+3) 

0031 

NUM=0 

0032 

IN'.JAP=0 

0033 

KN  =  1 

0034 

MULT-0 

.3335 

KM-0 

005b 

NRcC=0 

0O37 

NCCUNT-0 

0038 

C. 

0033 

C 

INITIALISE  I  COUNTER. 

0030 

ICOUHT=NI 1-1 

003 1 

c, 

0032 

c 

INCREMENT  I  COUNiER  AND  CHECK 

0033 

c 

IF  TOO  LARGE? 

0034 

1 1 

CONTINUE 

0035 

I  COUNT-- ICOUNT+1 

005b 

IF ( IC0UMT.aT.Nl2)  GO  TC  21 

0037 

c. 

0038 

c 

INITIALISE  J  COUNTER. 

0039 

JC0UNT=HJ1-1 

0100 

c 

0 1 6  ! 

c 

INCREMENT  J  COUNTER  AND  CHECK 

0132 

c 

IF  TOO  LARGE  ? 

0133 

12 

CONTINUE 

o  1  r,4 

JCOUNT-JSO'JNT+1 

0 ).  85 

IF ( JCOLNT. 3T. NJ2?  GC  T0  11 

P  1  :)6 

c 

0 1  37 

c 

INITIALISE  K  COUNTER. 

01-28 

KCCuNT=M.<  i  - ! 
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0139  C 
0110  C 
0111  C 
0112 
0113 
0114 
0115  C 
0116  C 
0117  C 
0118 
0119  C 
0120  C 
0121 
0122 
0123 
0124 
0125  C 
0126  C 
0127  C 
0128 
0129 
0130 
0131 
0 1 32 
0133 
0134  C 
0135  C 
0136  C 
0137 
0138  C 
0139  C 
0140 
0141 
0142  C 
0143  C 
0144 
0145 
0146 
0147  C 
0148  C 
0149 
0150 
0151  C 
0152  C 
0153  C 
0154 
0155 
0156 
0157 
0158 
0159 
0150 
0151 
0152 


INCREMENT  K  COUNTER  AND  CHECK 
IF  TOO  LARGE  ? 


13  CONTINUE 

KCOUNT=KCGUNT+l 
IF (KC0UNT.GT.NK2)  GO  T0  12 


HAVE  ALL  NUMBERS  READ  IN 
BEEN  STORED  ? 


IF(KN.LE.NUM)  GO  TO  17 


CALL  FFORM  FOR  NED  NUMBERS. 

14  CONTINUE 
NREC=NREC+1 

CALL  FFORMC IUNIT, IDENT.NREC, INEND,  INL*AR,NUM,  AA.  JDAT,  ININT) 
KN=0 


HAVE  ALL  NUMBERS  READ  IN 
BEEN  STORED  ? 


15  CONTINUE 
MULT  = 1 
KM=0 

16  CONTINUE 
KN=KN+1 

IF(KN.GT.NUM)  GO  TO  19 


IFCJDAT(KN) .NE.21  GO  TO  17 

SELECT  NEXT  NUMBER  AND  CHECK 

FOR  MULTIPLE  DATA. 

MULT= I NT ( A A ( KN ) +0 . 5 ) 

GO  TO  IS 

SET  UP  MULTIPLE  DATA  COUNT. 

END  OF  MULTIPLE  DATA  COUNT  ? 


17  CONTINUE 
KM-KM+1 

IF(KM.GT.M'JLT)  GO  TO  15 


NULL  DATUM  ? 


NCGUNT =NCO'JNT+ 1 

IF ( JDAT(KN) .EQ.0)  GO  TO  13 


STORE  NUMBER  IN  APPROPRIATE 
ELEMENT  OF  C  OP  IC  . 


L IM( 1 1 ) -ICOUNT 
L IM( JJ) “JCOUNT 
LIMCKK) =KCOUNT 
IFCiNlNT.EO.  1)  GO  TO  18 
C(LIM(n,LIM(2),LIMC3))=AA(KN) 
GO  TO  13  ^ 

18  CONTINUE 

IPUM=INT(AA(KN)+e.5) 

IF(AACKN) ._E.-0.51  IDUM=IDUM-1 
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0  1  33 

0154 

0155 

0156 

016? 

0158 

0168 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

017? 

0178 

0179 

0180 

0181 

0182 

0133 

0134 

0135 

0136 


IC(LIM(1),LIMC2),LIM(3)) = IDUM 
GO  TO  13 

C . 

C  CHECK  ENDING  INDICATOR  INEND 

C  AND  JUMP  TC  APPROPRIATE  LABEL. 

19  CONTINUE 

IFdNEND.GE.  IDIM)  GO  TO  21 
IF (NUM.NE.0)  GO  TO  20 
IF ( INEND. EQ.  1)  GO  TO  12 
IF( INEND. EQ. 21  GO  TO  11 

20  CONTINUE 

IF ( INEND.EQ. 1 . ANB.KCOUNT.NE.NK1)  GO  TO  12 
IF (JCOUNT.EQ.NJ 1 . AND.KCOUNT.EQ.NK1)  GO  TO  14 
IF (INEND.EQ. 2)  GO  TO  11 
GO  TO  14 

C . 

C  SET  ERROR  INDICATOR  I ERR  . 

21  CONTINUE 
IERR=NCOUNT 

IFUNUAR.NE.0)  IERR  —  IERR 


C  RETURN. 

RETURN 
END 

LIST  END  **** 
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0001 

0032 

0333 

0334 

(jo25 

033G 

0337 

C'038 

3052 


C 

C  FORTRAN  FREE-FCRMAT  DATA  INPUT  SCHEME 


5  :  FFGRM 
G  SUTLER  t.  J  PIKE 


MAY  '73 


C>WW<<^oKHo»T(oii))<'f<VM';  K^^^*****<^o!:>}:*.H<.KjMr^:iK}ko^st‘*)Ms>ioKHoK'+aK)M<^+-i<>K>,oMoK>WloWlo((*5i<HOk3(o|J 


-ftRAME'I’EftS 


0010 

C 

3311 

FFORM 

00  i2 

c 

nums a; 

0-313 

r 

£014 

*1 

1315 

n 

INPUT 

30  ib 

C 

8317 

c 

0O1E 

c 

IDF.NT 

0019 

c 

IN  1 1ST 

6020 

r 

INiiftR 

3921 

r 

TUN  IT 

0322 

c 

HR  EC  ■ 

3.923 

c 

0O24 

c 

3025 

c 

OJTr'U 

082S 

1? 

-:-o?v 

r 

ec  23 

c 

Rfi(4RI 

8229 

r* 

INEHD 

0-23.J 

INIJfiR 

0031 

r 

J  D  A  i  \* 

3032 

c 

tJo33 

C 

MUM  = 

033 -• 

r 

8L35 

r 

8-0  EE 

c 

IND  Id 

0O27 

c 

.  '’J'.  *'  fJ. 

r 

8079 

c 

INCC.-! 

c 

IN  Dm  7 

0041 

c 

0842 

c 

IMMG 

0643 

c 

INEXP 

0044 

c 

12! IT  ■ 

30  -C 

r 

cU  4L 

c 

INMUL 

00  a? 

c 

N-NL'U 

9340 

r 

INF- NT 

04  44 

u 

S3SJ 

CHORA: 

c~ 

— 

0,;C4 

C 

idL  ; 

1:021 

c 

CM  - 

80:  S 

t 

ID!  Ill 

U*  _l  ’ 

c 

I  Sit. 

L‘1j5  : 

c 

IDO  == 

j!.1  ’■  * 

r 

z * 

*  WARNING  COUNTER 
=  DATA  l.-.PUT  UNIT 


PARAMETERS 


13  HULL  SR  MULTIPLE 
]?.  r~  '•ftLUES  ECU  NT 


-  INDICATES  that  ftH  ALPHANUMERIC  COMMENT  HAS  STARTED 

-  IKE- ICrtTES  W-ETHSR  CURRENT  VALUE  IS  NORMAL, 

NULL  CP  MULTIPLE  5 A TP 

-  UT3ICACG  THAT  THE  CURRENT  CHARACTER  IS  ft  DIGIT 

-  iNDICfV  ES  C}|  EXPONENT  HAS  CHAFED 

-  INDICATES  TrAT  THE  NbfTiR  fSftC  c:  a’Mf-.PTOxS  ARE 
TO  HE  7  ■ITT-LISED 

-  IN’OiCH  -  I  HAT  ft  MULTIPLE  %'.5  HAS  STARTED 


inn: 


>  1  AST  DEE  I MEL  FG  ’riT  f  :-  JrP 


ELft'MC 

.'O',  S' 1ft 

I;  -  DISCS  0  TO  9 
'  DOLL  ftr. 

DOLOUR  PL-'-'SE 


ax 

< 
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70 :  * 

x  PL  -- 

-LOS 

e  Gsi 

L 

IP'!  = 

=0!-'T 

0'jSo 

c 

I3L  -- 

?LA;  rt 

2054 

c 

153  » 

3  IiNGLE 

0ii->5 

o 

1ST  - 

STAR 

»■:  2? 
CCT: 

r' 

L 

r 

OTHER 

-'HRT-IE 

L 

~  — - *  -  - - 

— 

eo?? 

C 

•  !  /’  l 

r  ■ 

in?  -  :CMTR0L 

ARAMETER  FOR  DARNING  RETURN 

‘■072 

r 

J  =  CO  OR  EOT  CH-. 

;\£C  T~  R 

w 

JCH<?0 '  »  CHAR? 

C7ES  STORAGE  ARRAY 

*'  l.i  /  ‘+ 

c 

JF  =  PTE  VIOL'S 

>VThC  ter 

ivfj/5 

2 

jpcscqj?  «  yp.s‘: 

INC  PC'S  IT  !0  s'  STORAGE  ARRAY 

0f>?o 

c 

jq  =. 

;-ffcsf.  jp 

U2?7 

c 

K  =  CCjNTER  F&* 

characters 

ttOf  P 

c 

KNuL  ■  POLITIC 

lt-  *  for  multiple  data 

»,  “;7‘  ‘ 

L 

KP  *  L3S  1  ISON-': 

LANK  CHARACTER 

b0:*"J 

o 

KP03  =  POSITHT- 

Zr  LAST  NON-BLANK  CHARACTER 

20;' 1 

r 

LE  =  ACCUMLLA  V 

R  PER  EXPONENT 

0032 

c 

L3 3  •-•  31  Or;  OF  EHT 

0033 

c 

LY  =  AXiMULftT; 

7.  FOR  CcCIM-L  PLACES 

RR?a 

c 

NCH  ’  0UV8ER  0.- 

CHARACTERS  PER  LINE 

0035 

c 

MUM  =  COUNTER  r 

OR  '-UMBERS  FOUND 

23  -£ 

r 

SG  »  SIGH  OF  HL 

■  r  -*■ 

f.  i? .?  V’ 

C 

y.  «  RCCUMULhTC; 

FOR  DIGIT?.  VALLES 

0039 

c 

00?? 

c 

0:;'J3 

c  #*#**#***• 

•;•!•  ^ajoSotoWR *t***r*v-  ■<  io^S^oKJtoKJMoialoioix#** 

9031 

c 

35322 

c 

0033 

SU6R0UTJNE  FFORMCL 

l-  IT,  IDEMT.NREC,  INEND,  INL'ftR . HUM,  AA,  JDAT,  INiNT) 

Or’3  4 

DIMENSION  JCN(Pl),  I 

D !  <  10),  JP03(SS) ,  fin  (40;' ,  jTAT  ’  40 1 

•1-35 

DATA  FDIC  1-/IH6/,  10 

I (2.'-*THl/,  ICK4)  'TH3/. 

GOj'o 

1  IDI'Si/TW,  I? 

1 ;  :61  /l.HO/,  IDI  (?)  /  :r$/»  iDI(9).'!HV, 

£«.  ;  < 

2  :5I  CPj/lPS/.  I1'1 

,  r- '! ,"  I  '  {'■)/ . 

0232 

3  ILL  H  /.  ICrf  ; 

“  iPT/«H./,  Il'II^'i .  JFL/1H+/, 

0035 

4  I E  X/ '  1  h z  S  p  I  r>  H  / 1  " 

V,  IDQ/'Ui"  ",  IDOL/ ’.AT/,  ISL/1H// 

c. . 

0  i  5 1 

r. 

REAP  LINE  OF  CHARACTERS 

0  1.02 

c 

S.  SET  ST-'ETIf'S  .LIES . 

8133 

TEN  *>10. 9 

0 1 3‘t 

NCH-S-j 

3135 

iC  !  M^l.NEM 

•O’.  1C 

JCH  i'M)  =  IKL 

2;  j :  ? 

JPC5CY;«1S. 

?  1  ?  3 

1  C*.'i- '  I-i’.'E 

P  i  39 

NC.*  P  1  -NCIvK 

0110 

JCMNCHPl;  *IB'_ 

oi ;  i 

RS'-I' .  :!JMii‘,590)  <JCH2rt),M-i,NCH) 

0112 

IHEliO-S 

0  1  1  i 

income 

01  i'l 

:nmll=c 

31 15 

i-\.jhs-=imdar 

5!  IS 

Oil? 

j*:ch 

H  I  !  3 


J.:  .! 


u;  Kf,e  J 

•1120  !>0 
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a  ip  i  ik it= i 

0122  C . 

S123  C  LABEL  1000  . 

3124  1009  CONTINUE 

3125  C . 

0126  C  INITIALISE  NEW  NUMBER 

3127  C  IF  FECUIRED. 

0123  IF ( IN  IT. ME. 1)  GO  TO  2 

3129  SG=  1.2 

0i 30  X=C.3 

3131  LY=9 

a  1  32  LZ-P 

3133  LE=3 

3174  L3G'-'. 

3 1 75  I  NT 21-1=0 

0136  INP;-->2 

0137  IN2-P-3 

0133  IHI'AT~  1 

0173  IN;T“0 

i\--Z  2  f.o  T'i-'E 

GUI  C . 

2142  C  SELECT  NEXT  CHARACTER. 

Si. ;3  K~y--i 

0144  .TO  JI 

0 1 45  ,TF  J 

0146  IF- J.NE.JBL)  KF  = J 

014?  ' F  .Mc.IOL)  KPQS -  1 


O  i  4C 

3 !  7,3 

c. . 

F-JL'iiLi 

6 1 50 

c 

IS  CHARACTER  FART  OF  COMMENT  ? 

0151 

IFaHC.CM.E3.li 

GO 

T!  5020 

0152 

c. . 

6 1 53 

r 

HAS  ENA'SC  CHARACTER  OCCURRED  ? 

0 1 54 

T'lNSND.GT.B) 

GO 

"0  5090 

0135 

0136 

r 

HAS  EV POMEM'  To"  rEI> 

5T5? 

if<:hev».gt.0) 

GO 

77  3009 

0150 

p 

0159 

c 

CHEF 4  r 7 "  MILL  TATA  AND  PRINT 

3150 

p 

L’P'3'  7  IF  Mot.Ti-wE  NJLL  DPTH 

0151 

r- 

u 

cp:lv  ef'ned. 

9152 

I  T  *n 

0 1  53 

IFfJ.rX'.ICM)  II 

=  1 

y  1 5  5 

iFI.i.r.L.  IV_)  n 

=  1 

G 1  5o 

7 ij.ftj.  ;iou  i 

’-1 

0  i  i:>6 

:F  il.-T.tv  GO 

'•fr 

w1 

0 1 5? 

'  F  ,;-5. 5 i  I 'Ml'  I 

:;r.C’ 

■-•0 

0  i  63 

:  ■  •.  II*  NJi'l.ti  j  i  ci) 

•"*  r 
l.j<  - 

O'  7 

1  •'  . 

:  r  v  Nflij..  .f3.il 

'jL 

1 

0179 

■.  ,20.  >  «0  5 

>  I 

•  t.-^n 

9.71 

T  5  _.r  '.1^05 

*  G 

.  1  -’r  7 

;a  ?'> 

!  r* 

0 1  73 

J’T-  "'.  HUT  !  1 

017  4 

M  •  :*  ’  '< 

0 1 7’5 

{  ~  r  :  ~  \  i  M_ 

C  :  7  b 

IrJ  T-  :0  J 

&  1  iV 

3  2  '7 :  ;■  - 

n  ■  ”ft 

c. , 

0  1  7- 

H‘>5  NULL  TATA  SEEN  FOUND  7 

/•  l  Jl* 

:  -  ■  ■  77". 33. 01 

GO 

o 

00 


<u 

< 
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0 1  . 

r 

01?2 

c 

CHECK  -OP  D I f5 i T . 

0133 

{MD’G=0 

0134 

DO  4  n=l,'.9 

0135 

IF (J .  EO.  ID !  CM) )  IMDIG'-l 

9136 

4 

CONTI. -SUE 

c 

0133 

r* 

KPS  NUMBER  STARTED  ? 

f'!36 

I“f  INM2JM.GT.02  GO  TO  2000 

6  1 30 

c 

•3  >.  3 1 

c 

IS  CHARACTER  A  DIGIT? 

IFC INDIG. EQ.0)  GO  TO  5993 

0L33 

r 

•3  1  4 

c 

START  NUMBER. 

Q’35 

;NKLT1=1 

s  J  1 

c 

013-/ 

c 

CHECK  UK ETHER  CHARACTER 

0  i  ?8 

c 

BEFORE  NUMBER  IS  STANDARD. 

'VJVf 

3rXJP.EC. IpL)  INNUn=2 

0330 

iFCJP.EQ.  7.10  IF  HUM- 2 

623. 

IrO.Jr  .EQ.  I'll)  SG---1.0 

J232 

IF(JP.SQ,  131.)  !NNU-W«CH 

0f>£*7 

IFCJP.EQ.  ICN;  ImHUM-NCH 

9:  u4 

1 1*  f„  JP .  £Q .  I  liQ ;  I  jSNUM*i'\Ch 

01.  3fj 

IFCJP.EQ.  T 30)  Ir^.-X-iCK 

0336 

IF ( JP.NS.  I-'T)  CO  TO  S 

C  ' :  .>  -• 

V-.  •  1 

2.--r:NT=  i 

92  93 

1.2  -  5 

0239 

LV--  j 

02  .0 

IFCJO.EQ.  IPL)  INNL‘U-3 

3*>  i  i 

IF  'J9.FQ.  ITT)  INNUi  !=3 

2  212 

iFcjo.ea. (Mi)  3G=-;.e> 

.0212 

IFC.’Q.EQ.  3  3L)  INNiJXCCK 

0214 

IFuIu.I2Q.ICM)  INNUi  :=NCH 

■3  "’IS 

:r ;jr. eg. idq>  innu:.=nch 

0216 

IFCjQ.S-J.  ISO)  iMNIJM-iSCh 

0  21? 

IF  I  ThU1. .  r.Q .  1 .  ANT .  JO . EO .  1ST) 

INNUM-2 

ti.;18 

5 

CONXMIJE 

f- "  1  Ci 

c 

0,723 

c 

CHECK  FOR  MULTIPLE  DATA. 

03':  1 

c 

°RIH7  '.JA'RN 1 33  IF  JON-COMPACT. 

0  : 22 

IF (  riMUL . MI . 1)  GC  70  6 

02  53 

KPCS=K-!NSJM 

J.;!  J 

1 3  C  Kr  ilJL .  EQ .  KPOS )  INNUrt»2 

U2>S 

.7  :ui.EQ.pp03-'  go  tq  s 

J3r TiSUMi-t 

1  .'.J 

02z:- 

t  ppT’r^fpJL 

0239 

GO  TO  399 

1"  23  , 

c 

0222 

PRINT  WARNING  FOR  NON-STANDARD 

:  2  j? 

r 

CUARACTEP  BE “OR 2  HUMBER. 

£.13* 

6 

CONTINUE 

0,233 

GO  TO  2000 

i7~^. 

GO  tq  ?oo 

6  2  jsti 

c 

02  ?5 

c 

LAEEL  2002. 

*1.240 

2009 

CONTINUE 
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9242 
0243 
£2'  44 
45 
024G 


BUILD  UP  NUMBER  VALUE . 


LY-LY+LL 
Ml  -X*  10.0 
DO  ?  M= 1 . 13 

iF(.J.ZO.  'V.(M))  X=FLOBTf.M-l)+Xl 


9-4? 

7 

CChTiNL't 

7j.l  iiD 

c 

IS  CHARACTER  AN  E? 

02  30 

C 

IF  SO  START  EXPONENT. 

•j15  ! 

IF'J.EQ. IPX) 

!NEXP=1 

cC52 

IF'J.EO. lEX) 

GO  TO  L  CUC 

c . 

0254 

r 

IS  CHARACTER  FIRST  DECIMAL 

0255 

C 

PC  I NT  OF  DIGIT  ? 

0255 

IFr  J.KQ. 

AND. INPUT.  ME.  1) 

IKJDIG=2 

0?n? 

IFCJ.EQ. IPT) 

INPUT- 

t\i  r.-- 

IFCINDIG.'EO. 

0)  GC  tq  40.00 

uj-jC4 

r.  .  . 

0250 

r 

IS  CHARACTER  AT  END  OF  LINE  ? 

0251 

IFCX.EQ.NCHPl)  GO  “0  4C30 

025? 

r 

32-:o 

2*- 

IF  CHARACTER  IS  DEC.  POINT. 

3234 

r- 

*w 

CHANGE  nijmoER  VALUE  ACCUMULATORS. 

0253 

IFfJ.EQ. IPT) 

!_/.-! 

0'2‘r  5 

IF (T. EC*.  f  T) 

LY=- 1 

02:' 

GO  TO  5000 

0?-  0 

c _ 

0259 

<•* 

u 

LABEL  3900  :  UPDATE  EXPONENT. 

0270 

3000 

CONTINUE 

02?  1 

02?  2 

r 

IS  TH IS  THE  FIRST  CHARACTER 

82?3 

i_ 

AFTER  E  ? 

32  ?4 

IF ( INEXP.GT. 

1)  GO  ~0  8 

P’VT 

r . . . . 

4  -  O 

c 

IS  CHARACTER  CLANK  +  OR  -  ? 

327? 

I !-  r  J .  E? . TC'  ) 

I HE  I  F -2 

O.L  l’ 

IFU.ZQ. 

P’EXP-2 

C-79 

If-  ( J . £Q.  I'  il> 

iNEXP --2 

3230 

IF  ( J .  PCI .  Ilil) 

!.SG=-  ' 

92".  1 

IF ( INEXP.GT. 

1)  GO  1 U  5000 

0242 

C. . . . 

0233 

0234 

0235 

C.?3(" 

025? 

82-0 

0233 

0230 

• 

0..32 
22  jji 
3234 
0235 
3V3S 
023? 

t*  'O'* 

0':?3 


BUILD  UP  EXPONENT  VALUE. 


3  COM  I  MU“ 

LEfl-Lv  10 
INDIGO. 

DO  9  M---1  5  0 

IF (J.NE. 'll! M) '  GO  70  3 
LE  - 1 1-  1  *lE  ! 

INENP-3 
IHlTG-t 
9  t  jNTK.JE 


„F !  TNE;'P.ER.  3-1  GO  T0  .0 
Ml? -2 
GO  TO  399 
10  CONTI  Ml'” 


IF  NO  DIGIT  HAS  OCCURRED  IN 
EXPCNEM.  PRINT  MORNING. 


<r 

o 

30 


<U 

< 

s 

H 
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13  CHARACTER  H  DIGIT  ? 

IF ( INDIG . EH .0)  GO  70  4S09 


0G33 
a  o'-!  4 

<j  :j 

0  ;  -  \ 

C . 

tv 

r 

IF (K.FQ.NC-iPl)  GO  T;j  4£G0 

GO  T'0  53O0 

IS  CHARACTER  A T  END  OF  LINE  ? 

'j  ■>  *}  c 

LABEL  4000. 

033:3 

4030 

oiTriUi 

\  ;  *  *  1 

..Til 

C 

CHECK  'Ji’ETHEH  CHARACTER  AT  END 

ostr 

r 

OF  NL'MEER  IS  STANDARD. 

7 

IF <J, S3.  ISO  GO  TO  !2 

5_»  ?« 1 4 

Ir  J.2~|.  !?_:•  GO  -r?  \-i 

<•  -  .  t~ 

u  i  o 

..  »■  vJ.L  -t  ■  L  •  1  i  )  <30  •  <v  -4 

L'  c  i  c 

IF-.J.’ii!.  IC'D  GO  T  .2 

*. 

1. 

0  3 1 V 

i ■:  3  .  t-.D ,  !  >< 3.i  GO  TO  ‘  0 

' 

•7)  5 1  j 

;rC,.cTi.  :  ,33  GO  r J  .2 

iiv  , 

T.2C .  13  jL)  gg  ~q  12 

^  Nvl 

V 

!1*.“  2 

.r  .7.  H2m  1" r;  (LG  1* 

b32i 

c' 

i?  CHARACTER  rS  *  PREPARE  FOR 

$ 

UT24 

c 

MULTIPLE  DATA. 

K 

93  20 

•FCINIHK..G3. !)  GO  'C  U 

033-6 

IF •:  IrVNr.GT.0-  GO  TP  a 

•33  .  7 

;f  ;  ;ne;t\gt.0.’  go  '•?  11 

IF("G._T.!;.3)  GO  Ti  11 

o.4 

ii*OI.LT,k!.5)  GG  TO  a 

o  j"  j 

llTfji.  ~2 

* 

0331 

lNDAT-2 

f 

3332 

r.TLvl«v. 

f'p  o  4  o 

i 

U  jj  r 

n 

LjM  Imue 

* 

33  35 

c _ 

j  3‘3 

r 

PRINT  HARM MS  FOR  NON-STANDARD 

-  -» v  • 

c 

CHARACTER  OT  Er'T  OF  NUMBER. 

5  *f‘" 

IF'K.EOC-I.A’-O.  I^G.ep,  p 

GO  TO  12 

V  ’  4? 

IMP-1 

0  5-0 

Ir  R"-K 

P3 

go  ’u  :ca 

’  '  J  c- 

*  n 

* »_ 

I.C'I-G  iNuE 

c . . . . 

*  J  V  4  - 

STORE  "  Li HER  '-A-LT  IN  ARRAY  AA 

0  .•  O 

r 

AND  SIT  AFR>'V  3D  AT  ”0  SHOLI  NULL 

c'  * 

p 

OS  i  .P_E  WV.* . 

03  i 

NdM-aUM-a 

.r-. 

LP-LbG  kLC-'-Y 

Tj  .•  1  *; 

r  A  •  f "JMi  '  f,  M’^TEN^-tP 

3  33 

TI  r-  T-l-'-Iun.  -IMDfif 

0.,.-:) 

T,-J? 

,*v:  •  / 

r. . .  . 

C  •<  50 

r 

ir  INTEGER  READ  ROUTINE,  CHECK 

3.  j-. 

r 

FOR  INTEGER  TALUS.  IF  NON- INTEGER 

;**’>  ^cr 

r 

FOUND .  FRIST  JARWING. 

* 

.  :0 

HTF  =  l 

_  j  •. 

•:.*•  unNT.sr.a  go  m  <? 

7  *  V£ 

.F'. :rT4T..,r.o.i  i;iF=o 

1 

0 .  ::  -< 

IF'  >ET>. v=. 0)  IMP  f 

G2  it- 

•  Fri.-i,'-. 3)  go  to  13 
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3351 

i  [  =  iNr..sa(suri)+0.5) 

23.-2 

! F  i (it  n .  LE .  -0 . 5 >  ! !  =  1 1 - 1 

0253 

IMP =3 

35: 5  i 

UG  TO  952 

l  00  J 

13 

CONTINUE 

07o6 

C 

•  •  « 

0552 

c 

SET  INITIALISATION  FLAG 

3358 

C 

FOR  ti’EU  Nl  TIBER . 

0353 

1 N I T  =  1 

3320 

c 

022  1 

r 

LABEL  5030  . 

0522 

50R0 

CONTINUE 

c 

032-1 

c 

CHECK  FOR  COMMENT  START  OR  END. 

0323 

IFCJ.SO.  1 7>Q 3  INCGM-- IHCCM+l 

3222- 

IF  (,T .  EQ  .152)  iNCCM- -NCQM+1 

032? 

IF<  INCOM.GT.  1)  5HCG!  1=3 

3323 

.F< INCOM.EQ. 1  -  Q0  70  14 

03?  3 

p 

0380 

c 

0231 

r-- 

lw 

CHECK  FOR  S  OP  /  AND  SET 

a;  32 

c 

ENDING  INDICATOR. 

3" ■  33 

IFCi.EP.  IS-_)  INEND "4 

0334 

IF ( INEND.  E:h4>  C-0  75  14 

0555 

IFM.NE.ID2U  GO  TO  14 

0235 

rjENT*  • 

9522 

IF (JF.NE. IDOL)  GC  TO  14 

9335 

T  iJ  :  -i 

?  ;ay 

IF(JQ.NE.U'GL)  GO  T2i  14 

0530 

I  MENU?. 

0331 

14 

CONTINUE 

0  vj2 

0533 

c 

IS  CHARACTER  AT  END  OF  LINE  ? 

0224 

IF(K.NE.NC-!Fi)  GO  ~0  1900 

0325 

c 

035(3 

c 

IF  multiple  data  definition 

033? 

c 

UNFINISHED,  PRINT  ulARNING. 

0333 

!F(TNffJL.Ni.l)  GO  73  12 

032-2. 

iiiF=; 

t>430 

J DPT (MUM) *  1 

042  1 

KF.-.T-KfflJL 

043? 

FI  TO  559 

3403 

15 

CONTINUE 

0i;.,'J4 

c 

0435 

c 

PRINT  FURTHER  L'AP.NING 

C436 

c 

INFORMATION  IF  RE  Hi  I RED. 

0i>3? 

if;:dent.f:3.0,'  go  "o  j? 

0435 

Ir'!  TDEN~.GT.999)  GO  TTJ  56 

0439 

ICUN  .irtP.E:-'.  INGHo)  GO  to  1? 

0*10 

16 

CONTINUE 

04!  1 

M'-hCriF  1 

r> *s .  o 
....  _  t 

21 

LHM 

o-’  ;3 

IF'JCH-.M)  .Et!.  IEL.PH1  .11. HE.  1)  GO 

TO  21 

0414 

t?-; IT? '  50..- -10)  <.  JChCil: . !•>!, fl) 

04  15 

!Nr Cn?I 

0416 

2? 

Iv-.'l- ' 

0-17 

: r i j pc?. ( n> . >2:2 . 1  al . - ?u .m.ne.d  go 

TO  21 

3413 

IF.  INuW.Nr.  il'lJcs)  iprc(50.543) 

(JPC'SOC1  ,  k  =  1 ,  M) 

0413 

!-•  iDENT._r.M-y)  G5  '-'0  U 

9  130 

<r 

o 

CO 


M 

x 

H 
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f\a? 

i ?  wu' 

-  it  V  ■ 

c 

04^o 

0-!-4 

c 

reiuf  : 

RETURN 

c 

e-’-2 : 

PRI!*7  l.sP£N  :  W> 

rF  REQUIRED. 

'  -i  “ 

49?  vj: 

i r c *  ,s .c;  i.T4 : 

i.  '50. 5/0) 

0429 

‘i-!73 

'  P  *.  1  1j  ‘  i  ;  .  :  * .  £  )  o .  < 

:o  1 

irEMT,M-::HCj:,  »VJ 

i-2 

IF  f  111'1.?'  .  >■>  •• 

!"■)  IDENT.NAEi  .  : .  •?.  j 

t.  ‘  1  j  c 

it" ;  lr'*’ .  r.  .-)  •JRi''? 

.  OO.  Ill  J  •  .  i  i 

’• 

[■  f. imp.  :r.  v'r.c 

•  I'T  -T.  NREO  ■  .F  u".  4.-.CX 

7?T-1I , JCnfKPKT) 

d -'35 

!? 

04/" 

».i:  '  TN;r  .-.7,  4:  ^=r- r~ 

04i? 

ji  t '  S  r'  ' ' 

0-,/H 

o  ':  -  ,  u.  ij. ;  ?  . : 

v .  5 .  >  isi.in? 

0435 

4/ 

c 

c 

500  fvt .'!,??( 30ft  1 . 

F0CMHT3. 

044  i 

•  _•  .  4*") 

5 1 C  C  ?  r^  !  1. 1  J  .  :  H  -  C 

r*  i  v 

/  a  w  , 

J'4'iw 

I  ("i  CHAR . 

,-o  smarts  sjmser  j 

0444 

523  FORMAT ( 111  , 10HFOR 

LPcNT®, 15, 6H» AFTER, 14, 8H  LINES 

&,  13, 

0445 

1  VH  CHAR . 

i2.i  IN  EXPONENT  ) 

6446 

530  FORMAT (  1m’  ,  10HFOR 

IDENT*,  15..  SH.  AFTER,  14,04  LINES 

4,  13, 

044? 

1  7H  CHAR .  ,2A1» 

12H  ENDS  NiiMiER  ) 

044Q 

54C  FORMAT < 1H  ,B0A1) 

044? 

H5U  f" .F:'..-  .  lr.  .c:r..'..d) 

:  156 

-Oi  r  .  -o” 

r ,  '5,  J  •  »  r.  I  J .  7 , 

i 

i  i •  *h  -oo'ir.  iiiTRi  :r\  r,«H  :-:  oo  ;ed 

1  -  .■ 

•  r.-  C:  FE-.4  WARM .HOC 

) 

.■wj3 

r  i: 

:  ->f. K 

LI 
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USERS'  GUIDE  TO  FORTRAN  FREE -FORMAT  INPUT  SCHEME 


The  scheme  enables  free-format  data  to  be  read  »into  various  arrangements  of 
the  prescribed  variables.  Each  routine  starts  on  a  new  line  of  data  and 

continues  until  either  the  specified  number  of  numbers  has  been  read  or  an  end  of 

read  symbol  (/  or  $)  is  found.  The  routines  are  designed  to  read  anything.  When 
ambiguous  or  unconventional  data  is  encountered,  the  most  reasonable  interpreta¬ 
tion  is  taken  and  a  warning  message  is  output.  A  read  failure  can  only  occur 
from  a  system  failure,  for  example  if  the  number  is  too  big  for  the  computer  or 

an  attempt  is  made  to  read  beyond  the  end  of  a  file. 

Various  calls  are  permitted  to  give  maximum  flexibility  in  the  choice  of 
input  variables.  The  scheme  is  programmed  in  standard  Fortran  IV  and  has  been 
tested  on  a  wide  range  of  computers.  Data  prepared  in  accordance  with  the 
latest  FOR77  recommendations  for  free-format  data  input  are  compatible  with  the 
read  routines. 

The  read  routine  calls 

For  real  numbers, 

n  real  numbers  READn 

(where  n  represents  an  integer  1  ^  n  <  9) 

one -dimensional  real  array  READA 

two-dimensional  real  array  READB 

three-dimensional  real  array  READC 

and  for  integer  numbers,  similarly, 

n  integer  numbers  READIn 

one-dimensional  integer  array  READIA 

two-dimensional  integer  array  READIB 

three-dimensional  integer  array  READIC 

The  parameters  are, 


READn 

(IUNIT, 

IDENT, 

IERR, 

XI 

9  •  •  • 

,Xn) 

READA 

(IUNIT, 

IDENT, 

IERR, 

A, 

11, 

12) 

READB 

(IUNIT, 

IDENT, 

IERR, 

B, 

11, 

12, 

NI, 

Jl, 

J2, 

NJ, 

IORD) 

READC 

(IUNIT, 

IDENT , 

IERR, 

c. 

11, 

12, 

NI, 

Jl, 

J2, 

NJ, 

Kl,  K2 , 

NK,  IORD) 

with  those  for  the  integer  reads  being  identical,  except  that  the  real  names 
XI, ...Xn,  A,  B,  C  are  replaced  by  the  integer  variables  Ll,...Ln.  IA,  IB  and  IC. 
In  detail  the  parameters  have  the  values: 
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IUNIT  input  device  number. 

IDENT  identification  number.  This  is  a  number  which  appears  in  warn¬ 

ing  messages  to  identify  which  read  call  has  generated  the 
warning.  It  is  also  used  to  control  the  type  of  warning  message 
output.  IDENT  =  0  suppresses  all  error  output  and  IDENT  >  1000 
provides  a  complete  list  of  the  symbols  read  by  the  read 
routine . 

IERR  is  a  variable  name  which  is  given  a  value  at  the  end  of  the 

read  call.  Its  value  is  set  to  sign  *N  where  N  is  the  number  of 
numbers  found  and  'sign*  is  +1,  except  when  a  warning  connected 
with  ambiguous  or  unconventional  data  is  given  when  it  is  set  to 
-1 .  It  should  be  noted  that  suppressing  the  output  with 
IDENT  =  0  does  not  affect  the  operation  of  IERR.  It  is  antici¬ 
pated  that  the  main  use  of  IERR  will  be  to  control  'batch' 
running  ( eg  to  prevent  programs  running  if  ambiguities  are 
present  in  the  initial  data  input). 

XI,..., Xn  n  real  variable  names. 

A,  II,  12  A  is  a  one-dimensional  real  array  name  and  the  numbers  are  read 

into  elements  A(I1)  to  A(I2) .  That  is  (12  —  11  +  1)  numbers  are 
read  into  consecutive  locations  in  the  array. 

B,  II,  12,  NI,  B  is  a  two-dimensional  array  of  dimensions  NI,  NJ  and  the 

Jl,  J2,  NJ  numbers  are  read  into  B(I,J)  where  II  ^1^12  and  J1  <  J  <  J2 

(see  IORD) . 

C,  II,  12,  NI,  C  is  a  three-dimensional  array  of  dimensions  NI,  NJ,  NK,  and 

Jl,  J2,  NJ,  KI,  the  numbers  are  read  into  C(I,J,K)  where  II  <  I  <12,  Jl  <J  < 

K2,  NK  J2  and  K!  <  K  «  K2  (see  IORD). 

IORD  controls  the  order  in  which  the  arrays  are  filled.  For  two- 

dimensional  arrays,  if  IORD  =  21,  the  subscript  I  varies  whilst 
J  remains  fixed  (that  is,  columns  are  read  in  matrix  terms). 

For  all  other  values  of  IORD,  J  varies  whilst  I  remains  fixed 
(rows  are  read  in  matrix  terms). 

For  three-dimensional  arrays  a  similar  system  is  used;  for 
example  IORD  =  132  means  that  J  (identified  with  2)  varies  most 
rapidly,  followed  by  K  (identified  with  3),  followed  by  1 
(identified  with  1).  The  other  significant  values  of  IORD  are 
as  follows: 
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IORD  = 


231 

implies 

JK.I 

order 

213 

implies 

JIK 

order 

312 

implies 

K.IJ 

order 

321 

impl ies 

KJI 

order 

values 

of  IORD, 

the 

order 

Form  of  free  format  numbers 

The  read  routines  will  read  numbers  in  any  of  the  following  forms: 


Sd,  S.d,  Sd.d,  Sd.  SdEd,  S.dEd,  Sd.dEd,  S.dEd 

where  d  is  a  single  digit  or  a  compact  string  of  digits,  S  is  +  ,  or  no 
character,  and  E  is  an  E,  E  blank,  E+  or  E-.  Numbers  are  normally  separated 
by  blank  or  comma  (delimiters).  Other  forms  of  number  of  delimiter  will  be  read 
and  accepted,  but  a  warning  message  is  given  so  that  the  user  can  check  that  the 
interpretation  taken  by  the  read  routine  is  that  intended.  Some  special 
characters  may  also  be  used  without  giving  warning  messages  as  is  detailed  below. 

The  following  additional  features  have  been  introduced  for  convenience  in 
using  the  system. 

(1)  Non-numeric  comments  can  be  inserted  in  the  data,  but  a  warning  will 
be  given  if  there  is  no  delimiter  between  comments  and  numbers,  eg 

1.23,  ABC,  2.45  will  give  no  warning 
but  1.23ABC2.45  will  give  two  warnings. 

(2)  Alphanumeric  comments  may  be  inserted  by  enclosing  them  between 
inverted  commas  (either  single  or  double).  In  this  case,  the  inverted  comma  is 
allowed  as  a  delimiter,  eg 

I .23"AB3"2.45  will  give  no  warning. 

(3)  Null  data  (£e  variables  retain  their  previous  values)  are  denoted 

by  commas  optionally  enclosing  blanks,  eg  123,,  ,  345  gives  four  numbers: 

123,  two  null  data  and  345.  ?!ote  that  comma  occurring  as  the  first  non-blank 
character  on  a  line  implies  a  null  datum,  but  that  comma  occurring  as  the  last 
non-blank  character  on  a  line  does  not  imply  an  extra  null-datum,  eg  ,  123,  345, 
gives  three  numbers:  1  null  datum,  123  and  345. 
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(4)  Multiple  data  are  denoted  by  i*  number  where  i  is  a  non-negative 
integer,  eg 

3*3.14159  means  three  values  of  3.14159 
3*  ,  means  three  null  data. 

When  *  is  used  in  this  way,  no  warning  is  generated. 

(5)  The  character  /  ends  the  read  call  forthwith,  and  any  further 
characters  beyond  the  slash  on  the  same  line  are  ignored.  Thus  /  can  be  used 
for  the  protection  of  further  data  or  as  a  method  of  retaining  the  existing 
values  for  all  further  variables  of  the  read  call. 

For  all  read  calls  except  the  two-dimensional  and  three-dimensional  array 
reads,  $  has  the  same  effect  as  /.  For  two-dimensional  array  reads  $  stops  the 
current  one-dimensional  sch-array  and  $$  stops  the  read.  For  three-dimensional 
array  reads  $  stops  the  current  one-dimensional  sub-array,  $$  stops  the  current 
two-dimensional  sub-array  and  $$$  ends  the  read.  Note  that  $  and  /  will  have  no 
effect  if  included  between  inverted  commas. 

(6)  Integer  locations  are  set  to  the  nearest  whole  number  when  supplied 
with  real  data  and  a  warning  is  given.  It  should  be  noted  that  integers  and 
reals  cannot  be  mixed  in  the  same  read.  Hence  a  Fortran 

READ  (IL’NIT,  FORMAT) K,  A(K),  B(K),  C(K) 

would  need  to  be  replaced  by 

CALL  READ4  (IUNIT,  IDENT ,  IERR,  XI,  X2 ,  X3,  X4) 

K.  =  LNT  (XI  +  0.5) 

A(K)  =  X2 
B(K)  =  X3 
C  (K)  =  X4  . 

Example : 

To  read  three  numbers  into  the  array  PTS  at  locations  4  and  6  from  input 
device  5,  use 

CALL  READA  (5,  1,  IERR,  PTS,  4,  6)  . 

If  the  characters  read  from  device  5  are 

X  =  3.14  Y  =  . 3E 1  Z  =  2 

then  a  warning  is  given  because  the  number  3.14  starts  with  an  =  symbol  (only 
space  or  ,  are  permitted).  The  form  of  the  warning  is 
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*****  DATA  read  warnings  ***** 


FOR  IDENT  =  1 ,  AFTER  1  LINE  AND  3  CHAR.  =3  STARTS  NUMBER 

X  =  3.14  Y  =  . 3E 1  Z  =  2 
+ 


PTS  (4)  to  PTS 


(.6)  are 


set  to  the  values  3.14, 


3  and  2  respectively. 
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